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Abstract 

We analyse supersymmetric contributions to Bg mixing and their impact on mixing-induced 
CP asymmetries, using the mass insertion approximation. We discuss in particular the 
correlation of SUSY effects in the CP asymmetries of Bg J/'4"P ^ind B^ 4>Ks and 
find that the mass insertions dominant in Bg mixing and B^ —>■ (pKs are and 
{.^2'i)LR,RL, respectively. We show that models with dominant (5f3)Li?,_RL can accomodate a 
negative value of S^Ksi agreement with the Belle measurement of that observable, but 
yield a Bg mixing phase too small to be observed. On the other hand, models with dominant 
(^23) LL,RR predict sizeable SUSY contributions to both AM^ and the mixing phase, but do 
not allow the asymmetry in Bd (pKs to become negative, except for small values of the 
average down squark mass, which, in turn, entail a value of AMg too large to be observed 
at the Tevatron and the LHC. We conclude that the observation of Bg mixing at hadron 
machines, together with the confirmation of a negative value of S^Ks 5 disfavours models with 
a single dominant mass insertion. 



1 Introduction 



The impressive performance of the B factory experiments BaBar and Belle provides the 
basis for scrutinizing tests of the standard model (SM) picture of flavour structure and CP 
violation in the quark sector, and opens the possibility to probe virtual effects from new 
physics at low energies. In the supersymmetric extension of the SM, a new source of flavour 
violation arises from the fact that, in general, the rotation that translates flavour eigenstates 
into mass eigenstates will not be the same for quark and squark fields, which implies the 
appearance of a new squark mixing matrix or, alternatively, that of off-diagonal squark 
mass terms in a basis where the quarks are mass-eigenstates and both quark and squark 
fields have undergone the same rotation - the so-called super-CKM basis. A convenient 
tool for studying the impact of this new source of fiavour violation is the mass-insertion 
approximation (MIA), which was first introduced in [1] and since then has been widely 
used as a largely model-independent tool for analysing and constraining SUSY effects in 
B physics. In the super-CKM basis the couplings of fermions and their SUSY partners to 
neutral gauginos arc flavour-diagonal and flavour-violating SUSY effects arc encoded in the 
nondiagonal entries of the sfermion mass matrix. The sfermion propagators arc expanded in 
a series in S — A^/m?, where are the off-diagonal entries and rhq is the average sfermion 
mass. We assume A^ <S m|, so that the first term in the expansion is sufficient, and also 
that the diagonal sfermion masses are nearly degenerate. 

Flavour-changing box and penguin processes as observed at the B factories are very 
sensitive to fiavour-violating effects beyond the SM, and the constraints on or measurement 
of nondiagonal squark masses will help to discriminate among various soft SUSY breaking 
mechanisms. In summer 2002, BaBar and Belle reported the first measurements of the 
mixing-induced CP asymmetry S^Ks iii (p^s, which at the quark level is 6 — > sss 

and thus a pure penguin process, which is expected to exhibit, in the SM, the same mixing- 
induced CP asymmetry as observed in B^ — > J/ipKs [2]. The experimental results, however, 
updated in summer 2003, paint a slightly different picture: 

Sj/^Ks = 0.736 ± 0.049 (BaBar & Belle) [3, 4] (1) 
5,^, - -0.39 ± 0.41 [5, 6] - I +0 45 ±0 43 ±0 07 g^g^, [g] 

Although the experimental situation in Bd — > (f)Ks is not yet conclusive, the deviation of 
S(j)Kg from Sj/^Ks n^^y constitute a first potential glimpse at physics beyond the SM, and it is 
both worthwile and timely to pursue any interpretion of these results in terms of new physics 
and to analyse their impact on future measurements to be performed at the B factories or 
at the Tevatron and the LHC, sec e.g. [9-12]. 

In the framework of MIA, the measurement of Sj/^Kc;, which is in agreement with the SM 
expectation, indicates that {6f^)AB, A,B = L,R, is small [13], whereas the result for S(f,Kg 
indicates a relatively large (^fs)^^- Furthermore, by including the constraints on (^23)^5 



from b S7, it was found [9] that, for average squark masses of order 500 GeV, only models 
with dominant (523)i'-R,Ki' can accomodate a negative value of S^Ks- 

insertions also determine the size of SUSY contributions to Bg mixing and, as a 
consequence, the mixing-induced CP asymmetries in tree-level dominated decays like e.g. 
Bs — > J/i^(t>-i which is one of the benchmark channels to be studied at hadron machines. 
Within the SM, the Bg mixing phase is very small, and consequently Sj/^^ expected to 
be of O(10~^). In SUSY, on the other hand, the third-to-second generation {b — > s) box 
diagram may carry a sizeable CP violating phase, which is described in terms of the same 
mass insertion (523)ab governing the CP asymmetry S^Ks- is therefore both important 
and instructive to analyse all 6 — > s transitions in the same framework, paying particular 
attention to the correlations between observables. This is the subject of this paper. 

Our paper is organised as follows: in Section 2, we recall the master formulas determining 
Bg mixing and the CP asymmetry in Bg J/i^<t> and discuss the SM expectations for the 
Bg mixing parameters and the experimental reach for Bg mixing at hadron colliders. In 
Section 3, we discuss the dominant SUSY contributions to Bg mixing in the framework of 
the mass insertion approximation. In Section 4, we present numerical results and discuss the 
correlation between the constraints from 6 57 and S^Ksi obtained previously in Ref. [9], 
and Bg mixing. Section 5 contains a summary and conclusions. 



2 Bg Mixing and the Mixing-Induced CP Asymmetry in Bg — > 
2.1 Master Formulas and New Physics Effects 

Let us begin by recalling ^ the master formulas for Bg mixing and the resulting mixing- 
induced asymmetry in Bg — > J/ip(f). Like for B^, the mixing angles p and q between the 
flavour and mass eigenstates in the Bg system can be expressed in terms of the B^ — B^ 
transition matrix element M12: 



q 



(3) 



where we have used AF^ ^ AM^ and AP^ <C P*"*. The resulting mass and width differences 
between mass eigenstates are given by 

AM, = -2Mi2, AP, = 2Pi2 cos Cs, (4) 

where = arg(Pi2/Mi2). P12 can be computed from diagrams with two insertions of the 
Ai? = 1 Hamiltonian and is dominated by the tree contribution. SUSY effects are very 



^ Here we use the convention \Bs)i = p|-B°) + and \Bs)2 = p\B^) — q\^g) where we define CP|P) 

-|P) and AMs = M2 - Mi and AF^ = Ti - 
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small, so to very good accuracy one can set 



12 



^ 12 ■ 



(5) 



In the SM, M12 is dominated by top quark exchange; the mixing phase is given by 



argMf^^ = 2arg(V,,V;:) = -2X'rj = 0(10-^). 



(6) 



In SUSY, there are new contributions to M12 induced by e.g. gluino and chargino box 
diagrams, which potentially carry a large phase and which we parametrise as 



Ml 



12 



which entails 



AM, = rli^Mf', Ar, ~ An^ cos 2/3„ 



(7) 



(8) 



assuming (5s ^ arg M^2 



rsM rpj^g above result implies that new physics contributions will always 
lead to a decrease of AF,, as was first discTisscd in Ref. [14]. 

Let us now discuss the effect of SUSY on the mixing-induced CP asymmetry in the 
tree-dominated decay Bg J /'4"Pi which is expected to be very small in the SM and hence 
highly susceptible to large or even moderate new CP violating phases. Although the final 
state J/il}(p is not a CP eigenstate, but a superposition of CP odd and even states which 
can be disentangled by an angular analysis of their decay products [15, 16], the advantage of 
that channel over the similar process Bg J/iprji^) is the comparatively clean, although still 
challenging reconstruction of the via — > K^K~ , whereas the r](^) is even more elusive. 
Once the CP- waves have been identified, the analysis of Bg — > J/'^(t> proceeds largely along 
the same lines as that of B^ J/ipKs, except for the fact that, in contrast to B^ mixing, the 
width difference APg cannot be neglected and entails a slight modification of the formula for 
the asymmetry. Without a separation of the final state CP-waves, the mixing asymmetry still 
depends on hadronic parameters describing the polarisation amplitudes ^o,||,± characteristic 
for the final state (^o,|| for CP-even and A±_ for CP-odd). One finds, assuming no direct 
CP-violation, 



sin AMd 



T(F: ^ J/^0) - r(i?o ^ J/iJcP) 



D Im 



j/^0) + r(i?o 



Podd 



-|- Im 



DFoddW+i^evenW 



sin AMgt 
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1 /AF, 
cosh 



+ Re 



■ , /AF, 

smh 1 

V 2 



(9) 
(10) 



and D encodes the polarisation amplitudes: 



as a hadronic quantity, comes with a certain theoretical uncertainty. Ref. [17], for in- 
stance, quotes D ^ 0.3 ± 0.2. 
The parameter p is defined as 

_ a(e: ^ j/^0) 

odd,even Mo\ 
Podd,even " ^^^o ^ J/V;0)„dd,even ^ ^ 

and can be computed from the AB — 1 effective Hamiltonian, yielding 

Podd,even ~ '^TF*TF~ ~ Codd,even (13) 
^cb ^cs 

with ^even = +1 and ^odd = ~1- Accordingly, we have 

~Podd,even — Codd,evene (14) 



2.2 Estimate of AM|^ and AVf 

In order to estimate AM|^\ one usually uses the ratio AM|'^'/AMJ^, in which all short- 
distance effects cancel: 

The remaining ratio of hadronic parameters has been calculated on the lattice yielding [18] 



BB^mi^f^ 



2 



:i.i5 ± o.o6^^o')^ 



where the asymmetric error is due to the effect of chiral logarithms in the quenched approxi- 
mation. In many SUSY models the dominant new contributions to Bd mixing involve transi- 
tions between the third and the first generation and are thus suppressed by the corresponding 
CKM matrix elements, so that B^ mixing is saturated by the SM contribution [11, 13, 19, 20] 
and we can assume AM^ — AMf^. AM^ is measured from the time-dependence of B^ 
mixing and is rather precisely known [21]: 

(AM<i)exp = (0.489 ± 0.008) ps'^ 

As for one has to use a value that is not contaminated by new physics. Stated 

differently, one needs a measurement of the angle a^^ or from pure SM processes. Various 
strategies for a clean determination of these angles have been proposed, see Ref. [22], and 
are expected to yield stringent constraints in the near future. For the time being, however, 
one has to resort to a different method and exploit the very basic fact that a triangle is 
completely determined by three parameters, which in our case arc the base, of length 1, 
the left side, which is determined by |V^6/K;b|, and the angle {3^^^ between the base and 
the right side. The essential assumptions that enter here are (i) that the determination of 




(a) (b) 

Figure 1: (a) Allowed region (shaded area) for the apex of the SM unitarity triangle, using 
the constraints from and sin 2/?. (b) AM^^ as function of 7^"^ as determined from 

(a). 



\Vcb\ and \ Vub\ from semileptonic decays is free of new physics, which is a model-independent 
assumption as these are tree-processes, and that (ii) /3 as measured from J/ipKs is 

actually - which, as mentioned above, is indeed the case in many SUSY models, but is 
a more mo del- dependent statement than (i). Using 

sin 2/3 = 0.736 ± 0.049 [3, 4] (16) 
I WKfel = 0.090 ±0.025 [21], (17) 

one obtains an allowed region for the position of the apex of the unitarity triangle which 
is shown as shaded area in Fig. 1(a). The allowed values of 7^^ are 45° < 7^^ < 100°. 
iVf^s/Vfdl can be read off the figure as a function of 7^^ from the right side of the triangle 
and translated into an allowed region for AMf^ as shown in Fig. 1(b), where we also include 
the error from Bs^fsJiBBaf's^- As can be seen from this figure, the current experimental 
bound AM5 > 13ps~^ [21] does not yet exclude any value of 7^^ between 45° and 100°. 

Let us now turn to AF^^. A recent estimate including NLO QCD corrections and lattice 
results for the hadronic parameters yields [23] 

—^ = (0.12 ±0.06). (18) 

s 

At present, there is no experimental bound. 
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Figure 2: (a) Correlation between Xs and 2/3^ for \R\ G {0.3,0.5,0.8,1,3,5} and argi? G 
[0, 27r], where R parametrises the new physics contributions to M12, Eqs. (19), (20). The 
numbers in the figure represent the values of \R\ and the circles and triangles indicate 
argi? = and vr, respectively. The value of argi? increases in the direction of the arrow. 
The perpendicular line is the current experimental lower bound of Xg- (b) New physics in 
Ar^. The numbers in the figure represent the value of \R\. lAF^I is always reduced by new 
physics and can even become zero. 

2.3 Observability of the 5^ - 'B^ Oscillation 

A convenient measure of the frequency of the oscillation is the parameter Xg-, defined as 

_ AM, 

Xg indicates the observability of the oscillation, which is governed by sin(xst/Ts); it is evident 
that the experimental resolution of rapid oscillations with ^ 1 is extremely difficult. The 
current experimental lower bound is Xg > 19; recent studies of the experimental reach of 
the BTeV [24] and the LHC [25] experiments indicate that Xg can be measured up to values 
Xg ^ 90 (note that the corresponding parameter in the Bd system, Xd has been measured to 
be 0.73). The performance of ATLAS, CMS and LHCb in analysing Bg J/'4"P has also 
been studied, which allows the determination of the correlation between the new physics 
mixing phase sin2/3s and the frequency Xg [25]. Although the sensitivity to sin2/?s gets 
worse as Xg increases, values of sin2(3g as small as (9(10~^) are within experimental reach 
for moderate Xg < 40. 

Let us now discuss the correlation between 2(3g and Xg in terms of contributions from 
beyond SM. For later convenience, we parametrise the new physics contributions as 



(19) 



which imphes 



2/?, = arg[l + i?], 



AM! 



\l + R\ 



(20) 



In Fig. 2(a) we plot the correlation between 2/3^ and Xg for different values of |-R| g{0.3, 

0. 5, 0.8, 1, 3, 5} varying the phase argi? between and 2tt. The value of AM^^ is chosen 
to be 25ps~^. The figure shows that the current experimental bound on Xg has already 
excluded some phase region for 0.5 < \R\ < 1. In view of the limitation of the experimental 
resolution, Xs < 90, it is clear that new physics can only be resolved if it is not too large, 

1. e. \R\ < 4. As for the mixing phase, 2/5^, small \R\ <C 1 will result in small 2/5^ that cannot 
be distinguished from the SM expectation, unless argi? is very close to zero or tt. For large 
SUSY contributions |-R| > 1, on the other hand, sin2/3s ~ 1 is very possible. 

Let us now discuss new physics effects on AF^. As discussed in [14,15], AF^ is always 
reduced by new physics due to the factor cos2/?s in Eq. (8). In Fig. 2(b), we plot AFg/AFg'^ 
in terms of argi? for different values of \R\. As can be seen from this figure, AF^ can even 
become zero for large values of \R\ and argi? = ±7r/2. 
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Figure 3: The time-dependent asymmetry of Bg 
given in the text. 
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J/ip(f> acc. to Eq. (9); parameters as 



Finally, let us discuss the effect of AF^ on the time-dependent asymmetry Eq. (9). In 
Fig. 3 we show the time- dependent asymmetry of Bg J / ipcf) for the parameter set AMg = 
25ps-\ AFfVFf* = 0.12, D = 0.33, \R\ = 1 and argR = n/2. Note that the maxima of 
the sin AM^t curve slowly decreases with t, which is the effect of the denominator of Eq. (9). 
Although this effect is rather small, it may be used to determine AF^ once experimental 
data become available in a sufficiently large range of t. 

3 SUSY Contributions to Bs Mixing 



The mass difference in the Bg system and the time- dependent asymmetry Sj/^^p depend 
essentially on M12 which can be computed from the effective AB = 2 Hamiltonian H^^"^. 



In supersymmctric theories H^^"^ is generated by the SM box diagrams with W exchange 
and box diagrams mediated by charged Higgs, neutrahno, gluino and chargino exchange. The 
Higgs contributions are suppressed by the quark masses and can be neglected. Neutrahno 
diagrams are also heavily suppressed compared to the gluino and chargino ones, due to the 
electroweak neutral couplings to fermion and sfermions. Thus, the B^-B^ transition matrix 
element is to good accuracy given by 



(21) 



where Mf^, Mfa and M^2 indicate the SM, gluino and chargino contributions, respectively. 
The SM contribution is known at NLO accuracy in QCD [26] and is given by 



Gp 
An 



{V,lVtsfSo{xt)ri2B[as{f^)]-'/'' 



(-^^bJ|,5i(//)) , (22) 



where SQ{xt) is given by 



Axt 



n.rf 



I 3xf In Xt 
A{l-XtY 2(l-xt)3 



(23) 



with Xt = {mt/mwY ■ Contributions from virtual u and c quarks are suppressed by the 
GIM mechanism. The short-distance QCD corrections are encoded in 7725 and J5, with 
7^25 = 0.551 and J5 = 1.627 [26]. 

Including gluino and chargino exchanges, H^^'^ takes the form 



ttAB=2 



(24) 



i=l 



where Ci{iJi), Ci{iJi), Qi{ij) and Qi{ij) are the Wilson-coefficients and effective operators, 
respectively, normalised at the scale with 



Qi 
Q2 
Q?, 
Q4 
Q5 



(25) 



The operators (5i,2,3 are obtained from Qi,2,3 by exchanging L R. 

In MIA, the gluino contributions to the Wilson-coefficients at the SUSY scale Ms are 
given by [27] 



Cf{Ms) 



at 



216m| 



24a;/6(a;) + 66/6(a;) (5^^) 



(26) 
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C'liMs) 
CliMs) 
CliMs) 
CliMs) 



a: 



216mi 



'216m? 



(27) 
(28) 



216m| 



2 { [504x/6(x) - 72/6(x)] (52'3)LL(52'3)i?i? ' l^'2Ux){5i^) U^i^) Rl) (29) 



216mi 



2 {[24a;/6(a;) + 120/6(a;)] (52'3)LL(52'3)Ri? - ^m^{x){5i^LR{5i^)RL]m 



where x = myw?^ and is the average down squark mass. Exphcit expressions for f'Q^x) 
and fQ{x) can be found in [27]. The Wilson-coefficients ^1,2,3 are obtained by interchanging 
L <-> in the mass insertions appearing in Ci^2,3- Note that the coefficient of the mass 
insertion (^23) ll{S23)rr in is much larger than the coefficients of the other mass insertions, 
which renders AMb^ and Sj/^^ very sensitive to these insertions. 

The chargino contributions to the relevant Wilson-coefficients, at leading order in MIA, 
next-to-leading order in the Wolfenstein parameter A and including the effects of a potentially 
light right-stop, are given by [19] 

= 7£l2E{l^^iri^.il1(^32)L + 2A(<53^)z.L(53"2)LL]L2(a;.,a;,) 



-2y,|yal%y;2[(53"2)ii((53"2)«i + A(53«)LL(53"jR^]i?2(x,,x,-,z) 

+Y,W^^V*V,,V;, [(5«2)L + m32)RL{Sl)RL\ ^X^, X,, z)] , 



Cf{Ms) = :r^J:U^2U,2VJ^Va {S^s2)Il + ^32) ll{S^i) ll Mx^Xj), 



12^1 



(31) 
(32) 



where Xi ~ rri^+fm^, z = mf /m? and the functions R2{x,y,z), R2{x,y,z), Lq^x^ij) and 
L2{x, y) are given in [19]. Ui^ and V^j are the unitary matrices that diagonalise the chargino 
mass matrix and Yt is the top Yukawa coupling (for more details, see [19]). Note that, 
neglecting the effect of the Yukawa couplings of the light quarks, the chargino contributions 
to C4 and C5 are negligible and that charginos do not contribute to C2{Ms) and C2{Ms) due 
to the colour structure of the diagrams; nonzero values at lower scales are however induced 
by QCD mixing effects. 

To obtain the Wilson-coefficients at the scale ~ m;, one has to solve the corresponding 
renormalisation group equations, which to LO accuracy was done in Ref. [13], with the result 



o ft 



(33) 



as{Ms)/as{lJ'). The coefficients b^'^' , cf"'"^' and arc given in Ref. [13]. 



In order to calculate M12, we also need the matrix elements of the effective operators 
Qi and Qi over Bs meson states. As usual, the matrix elements are expressed in terms of 



the decay constant /s^, using the vacuum insertion approximation; terms neglected in this 
approximation are included in a bag factor which is expected to be of order one. One has 

(EflQ^lB",) = -IruBjiBM, (34) 
- (^M^^h^-f'M (36) 

the matrix elements of Qi are the same as for Qi. The hadronic parameters Jb^ and 
Bi have been calculated on the lattice, yielding [28]^ Bi(mb) = 0.86(2) (^'^4), B2{mb) = 
0.83(2)(4), B3{mb) = 1.03(4)(9), Bi{mt) = 1.17(2)(t7), and B^iruh) = 1.94(3)(lf ); as we 
shall see in the next section, we do not need a numerical value for Jb^- 



4 Numerical Analysis and Discussion 

Let us now proceed to the numerical analysis of the impact of SUSY effects on AMb^ and 
sin 2/^5, which is most conveniently done by studying the ratio R, Eq. (19), of intrinsically 
supersymmetric to SM contributions to M12. We start with the gluino contributions, which, 
as discussed in the previous section, depend on the average down squark mass and on the 
ratio X — {mg/mq)^. In terms of the mass-insertion parameters ^23, it! can be written as 



R~9 = - aiirrig, x) [{5^^)Il + (^2^3)^^] + a2(m,-, x) [{5Is)Ir + {&rl 



d \2 



fd \2 



+a3{mg,x) (^23)^^(^23) RL 



a4[mq, X) 



{^2z)ll{52z)rR 



(39) 



with X — m^lrn^q. The coefficients ai{mq,x) depend implicitly on the Wilson-coefficients 
and matrix elements defined in the previous section. Let us pause here for a moment and 
consider what range of values for ^23 we actually do expect. Although our analysis is model- 
independent, we may nevertheless get some guidance for what to expect by looking at various 
SUSY models. For instance, with niq ~ nig ~ 500 GeV, the minimal supergravity model 
gives {5i^)LL ^ 0.009 + 0.001 i and rr,lr,rl ^ 0, while the SUSY SO(IO) model predicts 
(^23)rr — 0-^ + 0-^ ^ {^2'i)LL,LR,RL — [H]. Models with nonuniversal A-terms lead to 

^The overall sign is different from the one in [28], which is due to the different sign choice of the CP 
transformation; we chose CPIP*^) = +\P ). 
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((5^3)Li? - 0.002 + 0.005 i and {S^3)ll,rr,rl ^ [29]. We thus see that, although this is not 
expected to be true in general, a single mass insertion is dominant in many models. This 
implies that, for {623)11, rr {{623)lr,lr) dominated models, only the term proportional to 
{ci2{^q,x)) coutributes to R. We would also like to mention that (523)ab is already 
constrained by B{b — > 57), which yields 1(523)^^^,^?^! < ^ \{62s)lr,rl\ < 0(10^^) [30]. 

Numerical results for the x dependence of ai{mq, x) are given in Fig. 4, for two represen- 
tative values of the down squark mass, rriq = {300, 500} GeV. In order to obtain this result, 
we have set Ms = rriq and used the following input parameters: 

Vts = 0.0412, mt = (174 ± 5) GeV, a,(M^) = 0.119, 
mfe(mfe) = 4.2 GeV, fi = rrib, m,(2 GeV) = (100 ± 20) MeV. 

The impact of the theoretical uncertainties of mt and on Oj is very small, and also the 
variation with n ~ rub does not exceed a few percent. The main source of uncertainty of 
ai{mq,x) comes from the Bi parameters: although the factor Bi cancels in oi, the other 
carry a ~ 20% uncertainty from Bi/Bi. Note that Rg is independent of fs^- 

Let us continue with the discussion of the results depicted in Fig. 4. The solid and dashed 
lines refer to rriq = 500 GeV and 300 GeV, respectively. We see that all are monotonically 
decreasing functions in x and are by about a factor 3 larger for rriq = 300 GeV than for 
rriq = 500 GeV. Note also that ai{mq,x) becomes negative for large values of x. It is also 
evident that a4,{mg,x) is largest, in agreement with the remark in the previous section, so 
that the dominant contribution to Bs mixing through gluino exchange is expected to be due 
to LL and RR mass insertions. Although a2,3{mq,x) ~ (9(10) are also large, the constraint 
from B{b S7) on the helicity-flip mass insertions (^23)^^?,^^ renders their contributions to 
Bs mixing negligible. 
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As an explicit example for the relative size of the a^, we choose rriq = 500 GeV and x 
which yields 

R^(m^ = SOOGeV, x ^ 1) 1.44 [(^2^3)^ + {S',s)U + 27.57 \{6',,)l^ + {6^)1, 



- 44.76 



(Si)LR{Si)Rr] - 175.79 Us^) lM,) rr] . (40) 



Using the constraints from b S7, \{S22)lr(rl)\ < 10~^ and \{S23) ll{rr))\ < 1, it is evident 
that helicity-flipping mass insertions contribute C(10~^) to Rg, whereas single LL or RR 
mass insertions can yield 0{1) contributions. 

In Sec. 2, we have ah^cady discussed the dependence of AMs and sin 2/3^ on R, cf. Fig. 2(a). 
The constraint from b ^ S'~f implies that LR and RL mass insertions alone cannot generate a 
value of 2f3s larger than ~ O{10~^), which is too small to be observed at the Tevatron or the 
LHC LL and RR mass insertions, on the other hand, can result in sizeable - and measurable 
- values of the Bg mixing phase: for instance, {S23)ll = 1 x e*'^/'' yields AM^/AM^^' = 1.75 
and sin 2/5, = 0.82, while for {5^3)ll ^ (5|3)Ri? = 0.1 x e^^^/^o one finds AM,/AM|^ = 1.12 
and sin2/3s = —0.93. Note that for the same mass insertion, i.e. (523)ll = 1 x e'^/^, the 
smaller squark mass, rriq — 300 GeV, accompanied by a; = 1 gives about 3 times larger \R\, 
i.e. \R\ > 4, which is beyond the experimental reach at the LHC, as discussed in Sec. 2. 

Let us now turn to the chargino contributions. The chargino mediated processes depend 
on five relevant SUSY low energy parameters: rriq, mi^, M2, jJL and tan/3. With m^^ = 
150 GeV, niq = 200 GeV, M2 = = 300 GeV and tan (3 = 5, we find 

^ ^ 10-^(53"i)ll(53"2)ll + 2 X 10-^(^3-,)^ + 9.8 x l0-\5!i2) LLiSi;,) rl 

+2 X 10-\S^2)LLiS^2)RL + 2.4 X 10-\S^,)RLiS^2)RL + 5-4 X 10-\S^2)rl, (41) 

which is obviously much smaller than the gluino contribution. Even though the chargino 
contributions are very sensitive to the value of tan/3, an increase of tan/3 to 50 only entails 
an enhancement of the the first two terms in (41) from 10~^ to 10~^ - still not large enough 
to distinguish AM, and sin 2/3, from the SM prediction. 

Let us finally discuss the implication of the experimental data of the CP asymmetry in 
the Bfi (pKg process, S^Ks- As the underlying quark-level process is a 6 — * s transition, it is 
clear that this process is governed by the same mass insertions, (523)a_b- Since a possible hint 
of new physics may already have been seen in this mode, it is very interesting to analyse 
the implications of the experimental data on S^Ks for Eg mixing. Let us first recall the 
main result of the supersymmetric contributions to S^f^Ks previously obtained in Ref . [9] : the 
mixing CP asymmetry is given by 

sin 2/3 + 2R^ cos 5 sin(^^ + 2/3) + Rl sin(2^^ + 2/?) 



^4>Ks — 



1 + 2Rs cos 5 cos 9 a, + i?j 



(42) 



where 5 is the difference of the strong phase between SM and SUSY, but assumed to be 
5 = in the following (see [10] for a more detailed discussion). R^ is the absolute value of 
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rriq = rrig = 500 GeV 


Mass insertion sets 


Results 


5ll(rr) 


^RR(LL) 


Slr(rl) 


Srl{lr) 


AM, [ps-i] 


sin2/9s 




1 X e~*^/2 











10.7 





0.50 


1 X e-^'^/^ 











43.5 


-0.82 


0.59 








0.01 X e-^'^/^ 





24.9 





-0.36 








0.01 X e-^'^/^ 





25.0 


-2.8 X 10-3 


0.19 


1 X e-*'^/^ 


1 X 








4.39 X 10^ 





0.25 


0.1 X e-^^/^ 


0.1 X e-^'^/^ 








50 


0.87 


0.70 


= n);i = 300 GeV 


Mass insertion sets 


Results 


Sll{rr) 




Slr(rl) 


Srl{LR) 


AAf, [ps-i] 


sin2/9s 


S^Ks 


1 X e-'^/2 











87.6 





0.05 


1 X e"*^/^ 











115 


-0.98 


0.37 








0.01 X e-^'^/^ 





24.8 





-0.76 








0.01 X e-^'^/^ 





25.0 


-8.3 X 10-3 


-0.15 


1 X e-*^/2 


1 X e-*^/2 








1.26 X 10^ 





-0.52 


0.1 X e-^'^/^ 


0.1 X e-''^/^ 








128 


0.98 


0.65 



Table 1: Numerical results for AMb^, sm2(3s and S^pK^ for some representative values of 
((^32)ab {A, B = L,R) for rriq = m-g e {300, 500} GeV. 



the ratio between SM and SUSY decay amplitude and 9tj, is its phase, that is 

/ asusy\ 

R,e^'^ ^ [-jsirj (43) 

For rrig ~ rriq — 500 GeV, we obtain 

R^e'''^ ~ 0.23(5i J23 + 97.4(5^^)23 + 97.4(4^23 + 0.23(^^)23. (44) 

Considering the same constraint from 6 — > S7, we arrive at the conclusion that the LR or 
RL mass insertion gives the largest contribution to S^x^ while the LL or RR contribution 
is subdominant. In Ref. [9], we found that it is very difficult to get a negative S,j,Ks from 
LL or RR mass insertion dominated models without decreasing rUq. 

The most interesting result we would like to emphasize here is that mixing and S^Ks 
are dominated by different mass insertions: LL,RR and LR,RL, respectively. In Table 1, 
we present our results for AMb^, sm2(3s and Stpx^ for various sets of the mass insertions 
with rriq = rrig = {300GeV, 500GeV} ^. As we have mentioned above, the LL and RR mass 

^In this table, the phases are chosen to be negative so that S^Ks becomes less than Sj/^Xs (see the more 
detailed discussion in [9]). 
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insertions may lower the value of 5*^X8 and make it comparable to experiment if the SUSY 
masses are light enough. In this case, however, AM^ becomes so large that it cannot be 
resolved experimentally. On the other hand, although LR or RL dominated models can 
explain the experimental data of S^Ks and also predict AM^ ~ AMf"^, which is good news 
for the experimental side, in this case sin 2/3^ is too small to be observed. Thus, once the 
oscilation is seen with a large amplitude at the Tevatron or the LHC, all models with a single 
dominant mass insertion will he excluded. If the Bs oscillations are resolved experimentally 
with Xs < 90, the only surviving models predicting a negative SfpXs and an observable sin 2f3s 
and AMg, are SUSY models with combined mass insertions effects. An example of this class 
of models could result in, for instance, the following mass insertions {6^2) ab- 

|(54)«^|~0.5, 

\{5i)LR\ ^ |(fe)RL| ^ 0.005, 

arg[(523)i^i] - arg[(5^3)i?R] ~ 
arg[(523)Li?] - arg[(5^3)i?L] ~ 

which lead to: 

AM, ~ 40ps-\ 
sm2p, ~ 0.86, 
S^Ks ^ -0.7. 

Such nonuniversal soft SUSY breaking terms {LR and RL of order 10^'^ and large RR) are 
possible in models derived from string theory, as discussed in, for instance, Ref. [29]. 

5 Conclusions 

We have studied supersymmetric contributions to Bs mixing and the mixing-induced CP 
asymmetry of Bg — > J/ip(p in the mass insertion approximation, including constraints from 
other 6 — > s processes, in particular h ^ and B^ — > The SM predictions for these 

quantities are Sjj^^ ~ 10~^ and AM^ = (10 — 30)ps~^, depending on the value of 7. We 
have shown that in SUSY these predictions can change quite drastically, which is mainly due 
to gluino exchange contributions, whereas the chargino contributions to these processes are 
negligible. We find that values Sj/^^ ~ and AMg = (10 — 10^) ps~^ are quite possible. 
We also find that unlike their effects on the CP asymmetry of B^, — > (pKg, the mass insertions 
{^23)lr(rl) do not provide significant contributions to these processes, whereas {52z)ll{rk) 
imply a large AM, and sin2/3s. We have argued that a clean measurement of the S° — 5° 
oscillation and a significant deviation of S^k^ from Sj/^k^ would exclude SUSY models with 



TT 

"4' 
TT 

2' 
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a single dominant mass insertion, which predict either small oscillation and negative S^k^ 
or large oscillation and S^Ks — Sj/^Ks- 
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